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This  Interim  Technical  Progress  Report  covers  the  work  performed 
under  Contract  AF33(657)-1 1  265  from  October  1,  1963  to  January  1, 

1964.  It  is  published  for  technical  information  only  and  does  not 
necessarily  represent  the  recommendations,  conclusions  or  approval 
of  the  Air  Force. 

This  contract  with  Republic  Aviation  Corporation  of  Farmingdale, 

New  York,  was  initiated  under  ASD  Manufacturing  Technology  Laboratory 
Project  8- 1 1 1 ,  "Electro  Spark  Extruding. "  It  is  administered  under 
the  direction  of  Mr.  T.S.  Felker  of  the  Basic  Industry  Branch  MATB, 
Manufacturing  Technology  Laboratory,  Aeronautical  System  Division 
Wright-Patterson  Air  Force  Base,  Ohio. 

Mr.  J.H.  Wagner  of  the  Manufacturing  Research  Department, 

Republic  Aviation  Corporation  is  the  engineer  in  charge  of  the  project. 
Mr.  Gunther  Pfanner  is  cooperating  in  the  research. 

The  work  presented  in  this  report  was  done  by  B.  P.  Leftlieris  of 
Re-entry  Simulation  Laboratory,  Research  Division. 

The  primary  objective  of  the  Air  Force  Manufacturing  Methods  Program 
is  to  increase  producibility,  and  to  improve  the  quality  and  efficiency 
of  fabrication  of  aircraft,  missiles  and  components  thereof.  This 
report  is  being  disseminated  in  order  that  methods  and/or  equipment 
developed  may  be  used  throughout  industry,  thereby  reducing  costs  and 
giving  "MORE  AIR  FORCE  PER  DOLLAR." 

Your  comments  are  solicited  on  the  potential  utilization  of  the  infor¬ 
mation  contained  herein  as  applied  to  your  present  or  future  production 
programs.  Suggestions  concerning  additional  manufacturing  methods 
development  required  on  this  or  other  subjects  will  be  appreciated. 
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INTRODUCTION 


Capacitor  discharge  extrusion  is  a  metal  fabrication  process 
involving  application  of  the  effects  resulting  from  a  rapid  discharge 
of  capacitor  stored  electrical  energy  to  a  billet  in  order  to  produce 
work  in  the  form  of  an  extrusion.  The  electro-hydraulic  reduction 
of  ar.  extrudable  billet  rr  aterial  may  be  accomplished  in  a  number  of 
ways.  One  method  is  to  apDly  the  shock-pressure  energy  obtained 
from  the  capacitor  discharge  event  directly  Lo  a  billet  in  an  ultra-high 
pressure  fluid  environment.  Rapid  gap  discharge  of  electrical  energy 
in  an  appropriate  chamber  causes  a  shock  pressure  wave  to  strike 
the  billet  face  and,  upon  reflection,  to  create  a  dynamic  force  when 
added  to  the  latent  hydro-static  force  is  sufficient  to  effect  movement 
of  the  billet.  The  capacitor  discharges  are  intended  to  produce  a 
series  of  extremely  short  duration  intense  pressures  considerably 
in  excess  of  any  pressures  obtained  hydrostatically  in  mechanically 
available  in  tooling  of  comparable  size.  It  is  believed  that  the  short 
duration  of  the  pressure  pulses  will  allow  a  superposition  of  pressure 
without  the  corr  isponding  increase  in  the  design  strength  of  the  con¬ 
tainer  by  conventional  hydrostatic  criteria. 

With  the  billet  in  the  high  hydrostatic  pressure  medium,  an 
additional  phenomenon  of  the  capacitor  discharge  event  is  available 
for  extrusion  work:  namely,  the  electromagnetic  effects  derived  from 
the  rapid  discharge  of  electrical  energy  through  an  appropriately  con¬ 
structed  conductive  coil.  These  coil  discharge  rhet.omena  may  be 
applied  by  induction  to  the  face  of  a  piston-disc  interposed  between  the 
cou  and  the  billet  to  produce  a  shock  wave  in  the  high  hydrostatic 
pressure  environment  from  the  rapid  acceleration  and  sudden  stop 
of  the  disc  when  the  current  ceases  to  flow.  In  addition,  the  electro¬ 
magnetic  repulsion  may  be  applied  to  the  billet  face  directly  in  order 
to  realize  the  acceleration  forces  required  to  extrude  the  billet. 

A  further  method  of  capacitor  stored  energy  utilization  is  to 
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recover  the  energy  created  by  the  electro  hydraulic  discharge  and  to 
transfer  it  by  mechanical  means  to  the  billet  face.  Displacement  of 
the  billet  is  accomplished  mechanically  by  the  use  of  a  piston- ram  on 
the  billet  material  in  a  suitable  container  with  the  pressure  pulse 
energy  created  by  a  series  of  discharges  impinging  on  the  face  of  the 
piston  causing  it  to  move  forward  to  extrude  the  billet. 

The  advantages  offered  by  the  foregoing  approach  to  electro- 
hydraulic  extrusion  are;  (1)  elimination  of  billet  container  wall  friction, 
(2)  reduction  of  billet  die  surface  friction,  (3)  creation  of  pressures 
considerably  greater  than  hydrostatically  obtainable  allowing  extrusion 
ratios  beyond  those  obtained  hydrostatically,  (4)  reduced  container 
size  and  cost. 

Recognizing  the  attractive  advantages  of  the  application  of  this 
technology,  the  Aeronautical  Systems  Division  of  Wright  Patterson 
Air  Force  Base  has  awarded  Contract  No.  AF33(657}-11  265  to 
Republic  Aviation  Corporation  to  determine  the  production  potential 
for  extruding  steel  alloys  by  capacitor  discharge  energy.  This  will 
be  accomplished  in  the  course  of  the  program  by  conducting  experi¬ 
ments  with  capacitor  discharge  equipment  to  develop  the  most  suitable 
techniques  for  extruding  steel  sections.  The  program  consists  of 
two  phases  as  follows: 

Phase  I  -  Design  and  Extrusion  ol  the  Equipment 
Phase  II  -  Development  of  the  Extruding  Process 

The  methods  heretofore  described  will  be  investigated  tor  suit¬ 
ability  to  realize  *:.e  oojectives  of  the  program  of  extruding  a  2  inch 
diameter  billet  to  a  1/4  inch  round. 

This  second  quarterly  report  is  concerned  with  the  consequences 
of  the  explosive  effects  both  in  a  hydrostatic  medium  and  also  in 
pressure  vessel  walls  of  a  rapid  discharge  of  electrical  energy  across 
a  gap.  Analyses  have  been  prepared  and  original  equations  of  state 
derived  to  consider  the  effects  of  shock  pressure  wave  phenomena  as 
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they  will  be  applied  in  this  program. 

As  the  initial  experiments  place  the  billet  in  ar.  ultra-high 
pressure  fluid  environment  in  order  to  realize  the  advantages  derived 
from  the  hydrostatic  extrusion  process,  the  capacitor  stored  electrical 
discharges  will  be  used  to  furnish  the  additional  energy  required  to 
approach  the  target  extrusion  ratios  of  the  program.  For  this  reason 
examination  of  the  new  parameters  extant  becomes  vital  in  order  to 
understand  the  extrusion  process  at  the  billet  die  interface  to  predict 
ideally  to  what  degree  extrusion  may  occur;  and,  to  what  extent  the 
extrusion  process  can  logically  be  optimized  from  the  variables 
involved  in  the  electrohydraulic  approach.  This  work  is  presented 
in  the  first  section  of  the  report. 

Although  pressure  vessels  are  being  safely  used  at  hydrostatic 
pressures  to  300,  000  psi  using  well  known  standard  static  design 
criteria,  the  solution  of  the  problem  of  containing  similar  pressures 
in  a  vessel  with  a  superposed  shock-pressure  pulse  loading  was 
nowhere  readily  available  for  application  to  the  equipment  required 
to  attain  the  objectives  of  this  program.  For  this  reason,  an  original 
analytical  approach  was  pursued  that  has  resulted  in  elastic -plastic 
equations  of  state  that  take  into  consideration  the  effects  of  shock 
loading  in  ultra-high  hydrostatic  pressure  equipment.  These  derived 
criteria  have  been  applied  to  the  boundary  conditions  imposed  on  the 
container  design  by  the  electrohydraulic  extrusion  experiments  with 
tne  conclusion  that  a  container  design  that  will  accomplish  the  objectives 
of  the  program  is  feasible.  This  work  is  presented  ’r  rhe  second  section 
of  this  report. 
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NOMENCLATURE 
Part  1: 


p  =  density 

C  —  acoustic  speed 

P  =  pressure 

u  =  particle  velocity 

r  =  radius 

x  =  distance 

t  —  time 

0  =  angle 

F  —  frictional  force 

N  =  normal  force 

p  —  coefficient  of  friction 

m  =  mass 

=  axial  strain 
rQ  —  initial  radius 

xQ  =  initial  lengih 

v  =  Poisson's  ratio 

f  =  yield  stress 
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NOMENCLATURE  (cont'd.) 
Part  2: 


internal  pressure 
radial  stress 
time 
distance 

radius  during  the  wave  action 
particle  velocity 
hoop  stress 

density  during  the  wave  action 
density  before  the  arrival  of  the  wave 
radius  before  the  arrival  of  the  wave 
radial  strain 

hoop  strain 

da 

d« 

Poisson's  ratio 

yield  stress  from  uniaxial  tests 
Young'  j  modulus  of  elasticity 
plastic  modulus  in  the  bilinear  model 
shift  rate  (propagation  of  the  wave  front) 
acoustic  speed 

acoustic  speed  in  the  undisturbed  state 
propagation  velocity  of  the  wave  front  at  t= 
propagation  velocity  of  the  wave  at  rQ  =  a 


GENERAL  DISCUSSION 


From  previous  considerations  of  pulse  extrusion,  it  was  found 
that  a  pulse,  generated  by  the  discharge  of  a  capacitor  bank  at  one  end 
of  a  liquid  (preferably  water)  filled  cylinder,  a  pulse  of  215,  000  psi 
is  necessary  to  achieve  efficient  extrusion.  It  was  not,  however, 
proven  theoretically  that  such  extrusion  is  possible,  nor  that  a  cylin¬ 
drical  container  can  be  built  to  withstand  the  proposed  high  pressures. 

The  work  in  this  report,  therefore,  is  concerned  with  these  two 
problems  exclusively.  The  propagation  of  a  pressure  pulse  and  its 
ensuing  effects  require  a  different  theoretical  treatment  from  the 
existing  treatments  of  flow  (or  extrusion),  mainly  because  of  their 
short  duration:  in  gas  dynamics,  pressure  waves  represent  unsteady 
flow  as  compared  with  the  steady  flow  of  fluids  in  pipes  under  a  steady 
pressure  from  a  pump. 

The  formulation  of  unsteady  flow  dynamics  has  been  developed 
1  2 

extensively  ’  ,  It  is  not,  however,  a  simple  matter  to  apply  it 

because,  for  each  case  (with  a  given  geometry  and  a  given  medium), 
the  problem  must  be  formulated  anew. 

This  work  is  divided  into  the  following  two  parts: 

Part  1  Pulse  Transmission  and  Billet  Extrusion 
Part  2  High  Pressure  Transient  Loads  in  a 
Hollow  Cylinder 

Both  parts  are  generally  formulated  and  general  solutions  srt,  given 
and  discussed.  In  each  part,  however,  particular  solutions  are  given 
witn  assumptions  that  this  contractor  considered  to  be  reason:'  'o. 

The  results  are  sensible  and,  from  the  assumptions  made,  conservative. 
Nevertheless,  it  is  hoped  that  the  reader  will  form  his  own  opinions 
about  the  solutions  until  further  experimental  proof  is  provided. 
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PART  I  -  PULSE  TRANSMISSION  AND  BILLET  EXTRUSION 


The  general  case  for  pulse  extrusion  is  as  follows:  the  front 
part  of  the  cylindrical  billet  to  be  extruded  fits  into  the  tapered  section 
as  shown  in  the  sketch:  static  pressure  is  then  axially  applied  at  the 
larger  section  through  the  liquid  medium. 


The  capacitor  discharge  takes  place  on  the  left  end  of  the  cylinder 
and  propagates  to  the  right.  Upon  its  arrival  on  the  liquid-metal  inter¬ 
face,  it  is  partially  transmitted  and  partially  reflected. 

In  Part  I,  the  wave  interaction  at  the  interface  and  the  extrusion 
through  the  tapered  section  are  theoretically  analyzed. 
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SECTION  I  -  WAVE  INTERACTION  AT  INTERFACE 


Consider  two  cylindrical  bodies  of  the  same  diameter  which 
are  axially  in  firm  contact  with  one  another.  Consider  also  a  stress 
wave  travelling  along  the  axis  towards  the  interface.  It  is  required 
to  find  the  strength  ot  the  reflected  and  transmitted  waves. 


Cylinder  1 

//////////////. 

j 

|  \  — 1 

(J  Pi  Ct 

C? 

l 

1 

", -0 

Pj  =0,  p3 

r 

////////////// 

'///////////// 

Interface 


Two  condition:*  are  necessary  during  .reflection: 

(1)  The  velocities  of  the  wave  on  either  side  of  the 
interface  must  be  equal 

(2)  The  pressures  on  either  side  of  the  interface 
must  be  equal 

From  momentum  considerations  immediately  after  reflection,  we 
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obtain 


U) 


Pi  Cx  (U!  -U3)  sPj-Pi  *  Os  Cj  U3  -  Dj  Uj 


Hence 


P3  *  Da  C8  u,s  * 


-2  Pi  Cj  ui 
1  + -Pi  Cj 
Da  Ca 


It  is  obvious  that  when  pi-Ci  a  p3.C3,  Pa  *  Pi  (i.e.,  there  is  no 
reflection).  If,  on  the  other  hand,  P1.C1  <<:  p3.Ca,  then  P3  ■  2Pi 
(close  end  reflection). 


(2) 

(3) 


The  reflected  wave  is  given  by  (P3  -  Pi). 

<*-*>- ft 


Hence, 


(4) 


The  transmitted  wave  on  the  other  hand  is  given  by 
Pa  u3  C3 .  Substituting  from  eq.  (3)  we  have 
2ui 

1  ■  Pa  C» 

1+  Pi  6" 


u.  . 

a  p  •  c* + Pi  Ci 


(5) 


Thus,  if  the  amplitude  and  particle  velocity  of  the  incident  wave  are 
known,  the  amplitudes  and  particle  velocities  of  the  reflected  and 
transmitted  waves  can  be  calculated. 
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SECTION  II  -  EXTRUSION 


The  transmitted  wave  at  the  interface  propagates  through  the 
billet  towards  the  tapered  section.  It  is  now  assumed  that  the  die  is 
made  of  the  same  material  as  the  billet.  There  is,  therefore,  no 
reflection  as  the  wave  passes  through  the  tapered  section. 

Consider  an  element  of  thickness  dx  which  moves  under  the 
influence  of  the  passing  wave. 


From  the  momentum  equation,  we  have  LF  a  (6) 

where  ,  by  definition,  u  =  . 

N  sin  0  F  co3  0 

N  cos  0 

Also  F  =  Nfl  (?) 

EFi  (P  +  dP)  (r  +  dr)*ff  -  Pr*ff  -  2  (r  +  it  dr  (Np  cot  0  +  .  hi-  0) 

where  F 

The  mass  of  the  element  is  given  by 
dm  *  V  (r  +  dx  p 
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_  ,  ,  dr  ,a  dr  „ 

*  IT  (r  +  -T-)  p 

Thus  the  equation  of  momentum  can  now  be  written 

Zn  (r  +  4£-  }  p  »  (P  +  dP)(r  +  dr)  aff  -  PrBff  -  2{r  +  -^-)  n  dr 


(Np  cot  0  +  Nsin0) 


Transposing  and  neglecting  second  order  terms,  we  have: 


du_  m  2  _P_+  i_  dP 
ul  pr  p  dr 


-  2  cot  8  +  sin  0) 

r  " 


We  now  consider  mass  continuity  in  the  following  manner. 
Consider  an  element  initially  of  radius  ro  at  position  xo,  of  thickness 
dxo  and  density  p  o .  Under  the  influence  of  the  wave;  the  element  moves 
to  r,  at  position  x  and  its  thickness  changes  to  dx;  its  density  also 
changes  to  p. 

Thus  mass  continuity  is  expressed  as  follows: 

ff  ra0  dxo  p0  *  ffr*  dx  p  or  r*0  dr0  p0  a  rs  dr  p 


dx  a  dr  x  rBo  Po 

dxo  dr0  r*  p  (9) 

We  can  now  write  -jj-  »  since  we  considered  only  one  instant  of 
time,  and  "YY  since  we  considered  only  one  position.  Similarly, 

for  different  times.  Kence, 

aro  o  ro 


5u  2  _P_  ,  J_  jiP_  _N_  B 
?)t  Or  +  p  ?r  "pr 

wlic.e  Bap  cot  0  +  sin  0 
and  dx  3r  _  r8o  Po 

T5*  *  *5^  *  p 

Differentiating  (6)  with  respect  to  t,  we  obtain 


du 

dr0 


M 


(10) 


(11) 
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The  linear  strain  in  the  x  direction  can  be  written  as 


follows: 


Hence, 


ax  dxo  -  dxo 

£\  =  5xo _  =  'dx 

dxo 


Sxp 

ax 


dxo 


=  ft  +  1 


Combining  (6)  with  (8J,  we  obtain  the  following: 

|5-»|E  «{c1+  1)«  igl-fla- 
9xo  or0  1  '  r  p 


Hence, 


Thus 


4*1  =  -  (£c*  °°J  *(raP) 
^t  (rap)a  at 


^r0 


Mi 

at 


-t  o 


From  (9)  we  can  write 


1 

3r 


„  rap 


r0aPo  Sro 
Thus  (5}  can  now  be  written  as  follows: 


3u 

at 


rc  Vo  3r° 


-rr{p-B  »}■  ° 


(12) 

(13) 

(14) 

(15) 


(16) 


If  the  strain  in  the  two  directions  is  restrained,  the  strain  in  the 
ax’»!  direction  is  given  by: 

„  p  •>  p 
fi  *  E- 2v  r 


Thus 
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Substituting  in  eq.  (16),  we  ob':ain 


*U  r* _  _ E  .  2_fy  (2  -  B)  =  0  (17) 

"  r0*Po  "  (l-2v)  "  Sio  p  r 

Where  P  =  2fy  and  N  *  fy  substitutions  were  made. 


NOTE:  The  billet  is  compressed  continually  under  static 

pressure  so  that  upon  the  arrival  of  the  pulse,  the  material 
begins  to  flow  plastically.  The  total  available  energy, 

♦i _ fc ~c,  iz  giver,  by  z.)  the  static  fluid  prcsc-rr.  vvMcL, 

the  case  of  a  plasticized  solid  can  be  equal  to  twice  the  yield 
stress  and  b)  the  kinetic  energy  which  can  be  found  from 
the  momentum  equation  as  follows: 

From  the  momentum  equation  (A  P)  =  po  u  w 


where 


d£2 

dt 


dr 


3  U 


w  =  ii_=  41 

dt  dt  y  cTro ) 

dfrt 

if,  however,  »  o,  then  w  =  -g^- 
Hence  (AP)  =  u*p0 

Thus,  at  any  time,  (Ax)  A  AP  a  (Ax)  A  p0  v? 


where  A  *  Cross-sectional  area 
Ax  =  An  element  of  length 

Finally,  AP  _  us 
Po 

Equations  (15  and  (16)  form  a  system  of  partial  differential 
equations.  In  order  to  solve  them  simultaneously,  the  method  of 
characteristics  will  be  used.  The  two  general  equations  of  character¬ 
istics  of  two  dependent  and  '.wo  independent  variables  are  as  follows: 
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where  Ax »  Aj  t  Bi  ,  Bj  ,  Cx ,  Cg ,  Di ,  Dj ,  and  F i  ,  Fg  are  functions 
of  U,  V,  X  and  Y. 

Let  U  *  u 

V  *  fx 
X  =  r0 

Y  =  t 

Now  rewriting  (15)  and  (17)  in  characteristic  form,  we  obtain 


<\  «  "S  t  •  A  f  )  /  *V  «  \ 


0 


(&M&) 


r*  E 

royPo  IT 


L_  (*£l)  +  0  <2-B) 

2v)  \9r0J  °  \Bt  )  Por 


=  0 


where 


Ax  =  1 

Bx  =  0 
Cx  =  0 

Ox  =*-1 
El  3  0 


Ag  =  0 

Bg  =  1 

r  E 

~  r0a  Po  (l-2v) 

Dg  =  0 

Ea  m.Lh  (2  *  B) 

*  Por 


The  solution  of  the  system  is  given  as  follows: 

3 


4>  =  =  T 


dr  -  ^  r0^Po  (1  - 2v) 


r 

"ro*Po 


Po  (1  -2v)  '  V*p0 


E 

(l  -2v) 
E 

(l-2v) 


(18) 


.  ... !  r  \ 


1 


_  dfx  -  _2_  (2-B)  f  dt  =  0 

D0  (l-2v)  pr 


and 


du  Jr_\  E_  l _  df 

dt  \ro/^  Po  l-2v  dt 


dCx  _  2  fy  (2-B)  =  0 

P  r 


(19) 


The  partial  differential  equations  are  replaced,  therefore,  by  two  total 
differential  equations. 
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Now  let  p  =  p^  (constant  during  plastic  extrusion). 

Hence  +  1  a  ("r2-) 

We  shall,  furthermore,  assume  that  extrusion  takes  place 
in  the  plastic  region  with  E  *  0 

Equation  (15)  can  therefore  be  written 

du  _  2f„  (1-B)  a  0 

"3t"  "  Pc  r 

but  r  a  r  (r0,  t) 

Hence 

lr-(l rK-§fc)4? 


(20) 


(21) 


a  u  T 


r 

ro 


-L  JL 

po  l-2v 


a  U  1 


r  Jp.  '|Wv) 

»E.o.  -jf  -  If-  ■  • 

Thus  equation  (16)  becomes 


d3  r  2fv  (2-B)  1_ 

dts  Po  r 

A  familiar  and  exact  equation  can  be  derived  bv  rewriting 
eq.  (23)  as  follows. 


(22) 


(23) 


-§£-  d  a  HL-Bjjy  jdr_ 

dt  \dt  j  Po  y  r 

ru°_  /-ro 

I  du8  =  2  (2-B)fy  ( dr 

o  Po  Jrr 


Hence 


un8  _  4(2-_B)f  In  _ro_ 


Oo 


(24) 


15 


Equation  (24)  is  similar  to  the  steady-state  extrusion.  In  fact,  it 
is  identical  to  the  potential  energy  replaced  by  the  kinetic  energy.* 

Since  Ur,"1  is  a  positive  quantity,  we  deduce  that  0<B<1 
If  we  choose  B  »  0  (for  fX  a. 5,  8  *  39° },  then  we  have 

Up“  8fy 

P  0 


In  rp 
r 


(25) 


r 


n 

t 

„L 


-r 

i 


T 

i  i 

I 


*  The  mechanism  of  the  propagation  of  the  wave  is  as  follows: 
a)  the  wave  is  propagated  partly  with  an  elastic  steep  front  behind 
which  the  velocity  of  the  particles  is  Up,  and  b)  immediately 
afterwards,  the  decelerating  process  takes  place:  the  propagation 
velocity  of  the  deceleration  may  be  zero  if  £  ■  0  (i.e.,  the  flow 
is  plastic). 
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SECTION  II  -  EXAMPLE  OF 
PULSE  EXTRUSION 

Let  Pj  a  215,  000  and  fy  a  90,  000  psi  (fy  a  the 
yield  strength  of  annealed  4340  steel). 

Hence,  from  eq.  (5) 

2Pt 

U2  —  A  7*  x  P 

p*  •  r-ss  —  •» 

where  conditions  1  (subscript  l)  refer  to  water  and  conditions  2 
(subscript  2)  refer  to  the  billet  (steel  in  this  case). 

Ci  a  4800  ft/sec.  Pi  ■  62.4  lb/ft3 
Ca  a  1.8x  104  ft/sec.  pa  a  50C  lb/ft3 

Hence. 

Pi  Ci  +  Pa  Ca  *  (4.8  x  6.  24  x  104  +  1.8x5x10®)  -h  ^-c 


a  (30  x  104  +  9  x  10B)  2 


a  2.89x  10B 


Thus 

“a 


2  x  215,  000  4,3 
2.89  x  10®  2.89 


1 . 49  ft/ sec 


and 

P3  1  Pi  C,  a  2,8  x  10®  x  1.49  a  416,  000  psi 

These  are  the  conditions  immedi?tely  behind  the  wave  front: 
for  the  unloading,  Cs  a  us  a  uo,  and  Ps  a  AP. 

Hence,  from  equation  (25)  we  have: 

AP  In  _ro  and  with  fy  =  N  =  90,  000  psi 
8fy  '  r 

^2  i'qS  =  ln  *  This  equation  is  plotted  in  the  following  graph. 
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FIGURE  1 


SECTION  III  -  DISCUSSION 


The  result  of  the  mathematical  treatment  is  a  simple  equation 
similar  to  the  one  generally  used  l'or  static  extrusion.  The  two 
processes  are  not  directly  comparable,  however. 

In  the  static  case,  the  pressure  is  raised  until  extrusion  starts; 
this  is  the  pressure  necessary  for  extrusion  for  a  given  area  reduction. 

In  the  pulse  extrusion  case,  on  the  other  hand,  the  pulse  piopa- 
gates  through  the  billet,  causing  the  billet  to  move;  while  the  extrusion 
for  one  pulse  depends  on  the  pressure  amplitude,  the  final  extrusion 
depends  on  the  rate  the  pulses  are  generated. 

The  aforementioned  plot  of  AP  versus  ro/r,  therefore,  shows 
the  extrusion  for  one  pulse.  The  coefficient  of  friction  was  taken  at 
0.  5  and  the  angle  at  39*  (B>  Oat  0  ■  39*).  A  different  choice  of  values 
for  u  and  Scan,  of  course,  change  the  curve  shown  in  the  plot,  To 
find  the  rate  of  extrusion,  one  must  know  the  number  of  pulses  per 
second. 

n  q  p 

For  the  extrusion  of  ^s-y-  (i.e.,  »8  where  R0  is  the  radius 

of  the  billet  before  it  enters  the  tapered  section  and  R  is  the  radius  of 
the  extruded  section),  U  «  0.5,  0  ■  39*  and  AP  ■  400,  OOOpsi,  The 
extrusion  rate  can  be  found  as  follows: 


The  number  of  pulses  required 
to  extrude  tne  length  L  is 
given  by ; 

No.  of  pulses 

to  extrude  L  *  S®— /  ■£&— 

R  /  r 
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The  equivalent  length  of  extruded  stock  for  the  length  L  is  given  by: 


ffR*  li  * 


Hence, 


i,  .(|2)’ L  [;*•% 

■u^E-fr-)3  ♦  +•♦>] 

The  total  length  extruded  per  second  is  therefore  given  by: 


Ro 

tan  8 


[fSfiY  Ra  n  x  No.. of, 

JV_R  /  R  +  J  No.  of 


pul8es/sec 


pulses  to  extrude  L 


For 

-  8,  Ro  ■  1".  e  m  39\  ~~  ■  1 . 7  and 

^2  /_£ £  a  a  4.  7  pulses  to  extrude  1/ 

The  length  extruded  per  second  «  -g-j  x  +  ®  +  "i)x  ^^i_2LEHl££2Zl££) 

a  1  x  73  (No,  of  pulses/sec) 

.81  x4.7 

=  19.  2  x  (No.  of  pulses/sec) 

The  length  extruded  per  pulse,  for  ro  ■  1,  is  given  by: 

u  *  sS"8  EH**  -^+i] 

“  (2-9  +  1.7+1) 

a  5. 6 
.81 

a  6,95  in. /pulse 
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The  equivalent  length  of  extruded  stock  for  the  length  L,  is  given  by: 


'»•  11  ■  f  L  [‘  *  ft  * 
l,  -$)’  L  (**T5>  +(fe)il 

■a  \m'  ♦  *  d 


The  total  length  extruded  per  second  is  therefore  given  by: 

Ho  f/Re  \3  .  Ro  ,  T]  No.  of  pulaes/sec _ 

tan  0  |\  R  /  R  _J  X  No.  of  pulses  to  extrude  L, 


For 

4s2-  a  8,  R0  «  1",  0  «  39°,  —2—  »  1 . 7  and 

k  r 

S2  ^-f2  *  r~7  B  7  pulses  to  extrude  L' 

The  length  extruded  per  second  »  -g-j  x  (?*  +  ®  +  ^)x  pufses/sec) 

■  I  x  73  (Nc.  of  pulses/sec) 

.81  x4.7 

=  19.  2  x  (No.  of  pulses/sec) 

The  length  extruded  per  pulse,  for  r0  *  1,  is  given  by: 

*>  *  Stre  [(-?-)  ‘  *  r  *  ■] 

■  (2.9  +  1.7  +  1) 

*  5.  6 
.81 

*  6.95  in. /pulse 
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PART  2  -  HIGH  PRESSURE  TRANSIENT  LOADS 
IN  A  HOLLOW  CYLINDER 


The  stresses  and  strains  in  thick-walled  cylinders  under  internal  loads 
can  be  oaloulated  by  Lamp's  equations.  When  the  internal  load  is  transient, 
however,  equilibrium  of  foroes  does  not  exist;  Lame's  equation  therefore  ere 
no  longer  applicable.  Instead,  the  unbalanced  foroes  will  set  up  velocity  gradi¬ 
ents  in  part  of  the  oylinder  while  some  other  part  may  be  completely  at  rest. 

O 

Selberg  has  given  the  solution  to  suoh  problems  for  elastic  loads.  His 
results,  however,  are  not  explicitly  available;  it  is  not  possible  therefore  to 
either  modify  his  work  to  lnolude  plastic  loads,  or  to  use  his  results  directly 
for  plastic  deformation. 

In  this  report  the  general  problem  is  formulated  and  solutions  are  given 

4 

explicitly  for  both  the  elastic  and  elasto-plastic  cases.  Rinehart  hf.s  carried 
out  tests  with  thick-walled  cylinders  and  high  internal  loads.  His  results  indi¬ 
cate  agreement  with  the  prediction  given  in  this  report. 

The  solutions  can  be  used  for  the  analysis  of  stresses  in  shock  tubes, 
gas  guns,  and  pulse  extrusion  processes. 
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SECTION  I  -  PROBLEM  DESCRIPTION 


Consider  a  thick-walled  cylinder  of  inside  and  outside  radii  "a"  and  "b," 
respectively.  The  cylinder  is  filled  with  gas  or  liquid  to  a  static  pressure  P. 
Consider  also  a  pressure  disturbance,  generated  by  either  a  fast  accelerating 

piston  or  an  electrical  discharge,  travelling  from  one  end  to  the  other:  the 

da. 

disturbance  is  detonative  with  -jp  =  ®  at  its  front. 

The  pressure-time  profile  of  the  shock  wave  (detonative  disturbance)  is 
either  rectangular  or  triangular. 


In  the  triangular  case,  it  is  convenient  to  assume  an  equivalent  rectangular 
wave  with  the  areas  TPdt  and  JPdx  equal  to  the  corresponding  areas  under 
the  triangular  profile.  Thus,  the  momentum  and  energy  carried  by  the  wave 
are  the  same  for  either  profile.  The  triangular  wave  can  be  treated  separately, 
however,  without  difficulty. 

As  the  wave  travels  along  the  cylinder,  a  stress  wave  travels  radially. 
This  wave,  however,  is  not  cylindrical.  Its  velocity  of  propagation  has  two 
components:  v^^  the  velocity  of  propagation  in  the  medium  contained  by  the 
cylinder  and  v2,  the  velocity  of  propagation  of  the  elastic  wave  in  steel.  It 
is  assumed  in  this  analysis  that  the  angle  <t>  is  small  (i.e. ,  v2  »  v1). 
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The  problem,  therefore,  can  be  stated  as  follows:  Given  that  a  shock 
wave,  whose  magnitude  and  duration  are  known,  travels  along  the  inside  bore 
of  a  thick-walled  cylinder,  what  are  the  strain-time  relations  (both  radial  and 
hoop  strain)  in  the  wall  of  the  cylinder  if  the  yield  stress  of  the  mateiial  in 
the  uniaxial  test  is  known? 

The  resulting  longitudinal  wave  is  considered  out  of  phase  with  the  radial 
and  therefore  it  is  not  considered  in  the  analysis.  It  is,  however,  necessary 
to  assume  that  the  length  L  of  the  cylinder  is  considerably  greater  than  the 
inside  diameter  of  the  cylinder. 
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SECTION  n  -  PROBLEM  FORMULATION 


Consider  an  element  at  distance  r  from  the  center,  of  thickness  dr  and 
bounded  by  an  angle  d0. 


c  |  ■+  d<r, 


Thus,  when  the  wave  arrives,  the  element  is  accelerated  (or  decelerated) 
according  to  Newton's  third  law. 


where 


-dm.f-=SF 

2F  =  algebraic  summation  of  forces 

m  =  mass  of  the  element 

u  -•  particle  velocity 

t  =  time 
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From  Figure  1,  it  follows  tiwf: 


+  4  EF  =  -cr1Td0  +  (a1+dcr1)(r+dr)d6-2o^sin  y  •  dr 


but 

Hence: 


¥—»¥ 


2F  =  -al'r‘d0+a1-r*d0'<  a1*dr^i0+dal*i'd0+dP«dr'd0-c^*d0>dr 


W 


or  EF=  cr^dr«d0  +  da^r.dO  -  £^.d0*dr 

The  mass  of  the  element  is  given  by: 

dm  =  [Mfi±^±Mde]dr.p 

=  r.d0*dr.p 

where  p  is  the  mass  density.  Thus  the  momentum  equation  can  be  written  as 
follows: 

du 


or 


p  r*d0*dr*  =  o1*dr*d0  +do^*r«d0  -  a2*d0»dr 

du  °1  ,  ^1  a2 
"P  ~=~  + 


’  dt  "  r 


dr 


(2) 


but  only  a  specified  element  has  been  considered  so  far:  since  ^  will  change 

d  9  dq^ 

with  the  position  of  the  particle,  ^  =  gjf*  Since  is  for  one  instant  of  time 
dff1  3cr. 

only,  -jjr  =  "aT"  •  Hence,  equation  (2)  can  be  written  as  follows: 


The  continuity  equation  is  derived  from  the  fact  that  the  mass  of  the  element 
will  not  change  with  the  passing  of  the  wave.  Hence,  an  element  initially  lias 


*  The  equations  in  Part  n  are  numbered  independently  of  Part  I. 
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density  pQ  and  occupies  a  length  drQ:  during  the  passing  of  the  wave,  the  density 
changes  to  the  value  p  and  the  element  now  occupies  a  length  dr.  Thus 

dh  =  po-rQ.d0*dro  -  p-r*d0-dr  =  mass  per  unit  radius 


and 


3r__  Wo 
brn  "  p.r 


(4) 


Here  again,  the  total  derivatives  were  considered  for  one  element  at  one  instant 
of  time,  hence  they  are  now  expressed  in  partial  derivatives.  Th»  problem  can 
be  better  formulated  by  using  the  Lagranglan  representation  of  r=  r  (r0,t)  where 
rQ  is  the  initial  distance  of  the  center  of  gravity  of  the  element  from  the  center: 
thus ,  rQ  and  t  are  chosen  as  the  independent  variables. 


Equation  (4)  can  be  rewritten  as  follows: 
3r_ 


3 _ 3_ .  ^o 

3r  “  3r  3r  pt 
o  ^oo 

Substituting  (6)  into  (3),  we  get 


or _  3 


n 

p  at 


°ra2 


A  T*  3tri 

+  - i 

p  .r  3r 
*oo  o 


(5) 


Hence, 

au  r  *1  al"ff2  _  0 

(®) 

3t  p  •  r  3r  p.r 

nO  0  0 

Equation  (4)  is  now  differentiated  with  respect  to  time,  t,  and  hence, 

^•(p.r.|f)  =0 

0 

or 

^  „  du  ..  3r  do  ,  ^  dr 

p<r*  5F  +  r* aF  at  p‘  aiT  u" 0 

0  0  0 

(7) 

where 

u  =  |^  (by  definition) 

and 

p,r*  |r  +  rft  +  p,u  =  0 

(8) 
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Now  consider  the  radial  and  hoop  strains: 


-|£  dr  -  dr 

3ro  0  0  ar 

Radial  strain  €1  =  — jjj - =  Sr  "  1 


2jrr  -  2ffr 


Hoop  strain  «2  =  — 2~ — °  =  f-  -  l 


Hence  +  1  = 


£2  +  1  =  r" 
0 


Comhining  (4)  and  (9),  we  obtain 

<€i  +  l ,r  =  P0‘r o 

Differentiating  (11)  with  respect  to  t,  gives 
|-{<e1  +  l).p.r)  =  0 

(C1+i)p.|f  +  <€i+1)r-  |f  +  p*r.  ^L  =  0 


dp  p.r  J'l 

r*  a?-  -  P,u  “  (e1+l)  at 


.r  9el 


P'r*  3r  ■  (e^+l)  at 


Using  equation  (5), 


(D-r)2  3u  _  f[l_0 

p „.r  3r  €,+1  3t 
^0  0  0  1 


- 
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The  problem  is  thus  defined  with  two  partial  differential  equations: 


and 


r*U 

+  _L_  ^ 

St' 

P0*r0  aro 

dU 

-  _  A 

ar0 

aT  "° 

+ 


grg2 

p.r 


=  0 


(«> 

(12) 


For  the  solution  of  this  system,  rQ  and  t  are  the  independent  variables  and 
o^  and  u  are  the  dependent  variables.  Furthermore,  p=  p  (pQ,  o1),  r=r  (r0,  t), 
and  a2  =cr2  (r0,crj):  thus  the  coefficients  of  (6)  and  (12)  are  functions  of  a  com  ■ 
bination  of  u,  Oj,  rQ,  end  t.  Equation  (12)  has,however,  as  one  of  the  dependent 
variables.  This  indicates  that  an  additional  equation  is  required  to  solve  equations 
(6)  and  (12)  simultaneously.  The  equation  will  be  derived  from  the  equation  of 
state  of  the  material.  For  metals,  an  equation  of  stress-strain  will  be  sufficient. 


For  a  biaxial  stress  field  it  is  known  that 


€ 


1 


(N.B.:  The  cylinder  is  not  restrained  in  the  axial  direction.) 


(13) 


Whenever  the  stresses  extend  into  the  plastic  region,  the  problem  must 
be  divided  into  two  parts.  The  first  part  takes  place  in  the  elastic  region  as 
defined  by  Tresca's  criterion  of  failure,  i.e., 


assuming  that  the  stresses  cr2  an^  are  of  opposite  sign.  The  second  part  takes 
place  in  the  plastic  legion.  An  elastic-strain  hardening  bilinear  relation  for  the 
plastic  region  can  be  assumed  as  follows: 
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Hence 


and 


(ff-ffJ  =  En  (f-cj 


€  =!iE.u 

IP  Ep  MP  Ep 


In  general  we  can  write 


3c,  j  3a^  „  3crg 

•5t  =  e  at* "  e  "sr 


But 


12= 

E  ro 


1  5g2  1  3r  _  u. 
E  3t  ”  rQ  5t  “  r0 


Hence 


3t 


ro 


Equation  (18)  is  now  substituted  in  equation  (12),  giving 

3o, 


M  _  i._i  +  JL  u  =  0 

3r  E  3t  r 
0  o 


Rewriting  equation  (6)  we  also  have 

_L_  Oh  +  °JL?1  =0 
3t  P0*r0  3ro  P'T 


(14) 

(16) 


(16) 


(17) 


(18) 


(19) 
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The  solution  of  equations  (19)  can  be  afforded  by  the  theory  of  characteristics 


The  general  equations  are  as  follows: 


where  the  coefficients  A^,  B^,  C^,  D^,  Ej,  A„,  Bg.  C^,  D2,  E2  are  functions 
of  U,  V,  X  and  Y.  Now  let  u  =  U,  =  V,  X  =  rQ  and  Y  =  t.  Thus,  rewriting 
equations  (19),  we  have 


The  coefficients,  therefore, are  as  follows: 


Ai  *  1 

a2  =  ° 

Bi-0 

B2  =  1 

O 

II 

H 

C2  =  - 

“i--F 

J3 

II 

O 

II 

CM 

W 

The  solution  is  now  given. 


Povo 

pr 


(20) 
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and 

du  +  ♦  ~ —  dcr.  +  [M*  +  ~  ”  ♦Idr  *  0 

P0*r0  1  Lr0  p.r  *Jaro 

(21) 

where 

dt 

(22) 

o 


The  method  of  characteristics  has  thus  replaced  the  two  partial  differential 
equations  with  two  total  differential  equations.  Equation  (20)  has  the  units  of 
velocity  and,  In  fact,  is  the  propagation  velocity  in  the  Jagrangian  designation. 
Equation  (21)  shows  the  condition  of  the  metal  behind  the  propagating  front. 
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SECTION  m  -  BOUNDARY  CONDITION 


To  find  the  boundary  condition,  or  the  motion  of  the  interlace,  at  time 
t  =  0,  consider  an  element  dr  as  before.(See  Figure  1). 

Hence  dm  =  A  F 

and  dm  =  p.r*d0.dr  =  prQdfldro 


♦  +  ve  AF  =  -  c^-dS-r  +  (o^  +do1)  .(f-«ir)*d0-2cr28in  ^ dr 
=  +dtr1.r*d0  +a1»dr.d0  -  Og.dO’dr 


Therefore 

p.r.dr  =  da^r  t^-dr  -  or2.dr 

or 

dr 

p.r  g^  du  =  dcr^. r  +(ax-CT2)dr 

since 

r  =  r  (r0.  t) 

Thus 

5T=u+(l^)C(from^=|f J  =  u+(i+()C) 

0 

or 

! 

p.rud'i  +  p.r  (1+  e)Cdu  =  dcr^r  +  (cr^-Og)  dr 

Now  let 

dr-*0,  i.e.,  r- rQ  (dr  =  d  (r- rQ) ) 

Hence 

€  -  0  and  u  -  uQ 

p.u0duo  +  p.C0  duQ  =  do1 

p(u0  +  CQ)  du  =dcr1 

but 

uQ  +  CQ  =  wQ  (propagation  velocity  at  t  =  t0  ) 
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Pwoduo  =  dcrl 
and  for  a  sudden  change, 

p.Wj.u^  p  (23) 
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SECTION  IV  -  ENERGY  CONSIDERATION 


The  conservation  of  energy  in  the  cylindrical  stress  wave  front  can  be 
written  as  follows: 

u12p1«rr1dr-p|fra.dru2a+Bl0-Beg 

It  can  be  assumed  that  Ej08SCS  (strain  energy  losses)  is  negligible  for  moderate 
strains.  Thus  it  can  be  stated  that  for  all  stress  wave  cases  where  excessive 
plastic  "stretching"  does  not  occur,  the  kinetic  energy  per  unit  radius  remains 
constant.  Hence 


2  2  dE 

U1  rl  =  u2  r2  <i,e“  df  =  constant) 


(The  density  iu  this  context  is  assumed  constant.) 

It  is  thus  possible  to  find  the  velocity  at  any  radius  once  the  velocity  of 
the  bare  surface  (boundary)  is  known.  From  equation  (23)  u1  is  found  as 
follows: 

w1  =-  Cj  (the  acoustic  speed)  since  C^^  »  u1 
P  =  internal  pressure. 
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SECTION  V  -  SOLUTION  O?  EQUATION  21 


A.  GENERAL  SOLUTION 


It  was  found  previously  in  equation  (20)  that: 


and 


^2 

dt 


du  +  ♦  — —  - 

P0ro 


ala2 

p.r 


(21) 


These  equations  can  be  rewritten  as  follows: 


or 


du  .  1 

/  E  .  r  .  ^1 

dt  E  « 

7  p  r  dt 

du  ,  C 

^1  .  u.u.c 

dt  E 

dt  ro 

"o  ‘0 


o 


Before  any  attempt  to  integrate  the  above  equation  is  made,  it  is  important  to 
list  the  events  in  sequence. 


•  At  t  -  0,  the  internal  pressure  P  is  applied  suddenly. 


From  equation  (23).  u,  is  found.  The  velocity  and  density  changes  immediately 
behind  the  fiv.nl  are  sudden  and  take  place  almost  instantaneously  upon  arrival 
of  the  shock  front.  This  is  a  characteristic  of  «hock  waves  generally.  jU'e 
other  hand,  the  hoop  and  axial  strains  will  only  take  place  with  time,  as  the 
individual  particles  start  moving  under  velocity  u.  Thus,  from  equation  (U) 

•p-  =  1  (at  t=0),  the  radial  strain 


(v 


;  easily  seen  that  while 


is 


only  a  function  of  density  changes 


*o 

P 


Furthermore,  the  sudden 
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radial  change  will  be  restrained  both  in  the  hoop  and  axial  directions. 
Hence,  at  t  =  0, 


.2  gl 


.2  gl 


C1  =  E  '  ^  E  "  ^  E 
The  elastic  loading  will  thus  be  given  by 

V°2  gl-^gl 


or  ox  (l->t)=ay  (25) 

(i,  e. ,  c3  =  «2  =  -  fi  ^ .  Since  this  strain  is  restrained,  it  affects  the  radial 


gl 

strain  as  (H  jr)) 

and  el  =  Y  (1  “  2**2) 
Thus,  ~  =  1+y  (1-2.U2  ) 


(26) 

(27) 


at  t  =  0  and 


f  is  given  by  the  following  equation.* 


Vu 


(28) 


*  The  continuity  equation  across  the  wave  front  can  be  written  for  t  -  0  as 
follows: 

p  (C0-u)  -  p0  CQ. 

P0  „  Co-u 

Hence  j=C*  -g— 

Thus,  if  CQ  and  u  behind  the  wave  are  known,  £  can  be  calculated. 
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Immediately  behind  the  loading  wave  (t  =  0),  however,  a  triaxial  stress 
field  exists  and  the  maximum  radial  strain  that  can  be  applied  elastically  is 


al  ^  l-/i 


The  axial  and  hoop  stresses  are  compressive  at  t  =  0.  For  ct. 


>3l 

1  l-n  ’ 


the  ratio 


P0 

—  is  constant,  i.e.,  plastic  deformation  takes  place  at  constant  density.  At 

P 

t  =  0,  — -  =  1.  Equation  (29)  shows  that  the  elastic  limit  at  t  =  0  and  ~  -  1 
o  O 

is  higher  than  for  a  static  case  wherein  the  application  of  is  associated  with 

corresponding  hoop  stress  (of  opposite  sign).  Thus  the  plastic  deformation  be¬ 
gins  atffj-ffg  =  oy  . 
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For  a  compressive  internal  pressure,  is  negative  and  Og  is  positive. 


Thus 


-  KM<t2I=  ky|  < 


(30) 


/shock 

With  these  considerations  in  mind,  equation  21  can  be  written 

***?“*& ‘GO-*  <su 

It  is  convenient  to  express  the  different  variables  in  dimensionless  form  as 
follows: 

s-j 

0  0  0 


-  1L  s  _  E_  r  . 

ffl "  En  E  “  E  V 

O  0 


Hence,  by  substituting  in  (31), 


—  +M-O.C  +(51-52)C  (7)  =  0 

-  ITT  -  U~  f«L  £.  =  c  /e  — 

V  p_  r„  V  E„  p„  r  oV  r„ 


dt 

C 


^o  o 


Since 


r  =  f(r0,t) 


"o  Po  'o 


and 


dt 


=  5+CfC 


v 

Now,  consider  the  case  where  Oj  =  . 


(32) 


(33) 
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The  "stretching  out"  at  t>o  takes  place  in  the  plastic  region.  Since  C  =  C  =  0 
in  this  region,  equations  (32)  and  (33)  can  be  rewritten  as  follows: 


du 

dt 


+  <V&2  >*(f)  =0 


(34) 


B.  LOADING  WAVE  SOLUTION 

Let  f-  =  1  +  £ 
o 

where  £  «  1 
dt 

„  dV  du 

dt  dt 

d2/ 

Therefore,  — 5-  +  .rry.  =  0 
dt^  \1+Sl 

where,  from  Tresca's  criterion,  5  =  (a,-  a„). 

y  1  * 

~f'dy  C/  =  -5yC.  (35) 

The  quantity  dy  is  negative  since, we  are  considering  compression.  T^ns  (35) 
can  be  rewritten 

,Z£ 

^f+ayC,r=9yC  (36) 

The  solution  of  this  equation  is: 

£  =CX  eut  +  C2 e_a)t  -  1  (37) 
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40 


where  E  =  -  P*88t-—  .  Equation  (33)  can  therefore  be  rewritten  as  follows: 
^elastic 


dt 


o  rD  -  V  Fl/2 

P-a  +  £_£.c  =  u+£i.=  u  +  C  — - 


r_\l/2 


Xh) 


(41) 


Let  ~  =  1  +  £ ,  where  £ «  1.  Hence 
o 

—  =  u  +  C  E  1/2  (1  -  1/2  £) 
dt 

Equation  (32)  can  also  be  rewritten  as  follows: 


(42) 


+1/2C  El/2  + ME 1/2  (1+1/2  £)(—  +  1/2CE1/2^ 
dt  at  'dt 

-CE  1/2)  +  FC(1-  t)  =  0 
Rearranging  equation  (43),  we  have 

lif +(  1/2 C+M)  E1/2  |f-ayC^=C(ME-5y) 

Equation  (44)  has  the  following  general  solution: 

‘  2  +  2  J*‘ 


(43) 


(44) 


£-c,. 


■2+4o„C 


_S .I  IrT^.r. 


'y  >  O  O 

y  +  C2e  2  2  ' 


4py'  (uE-M 

\  ^  > 


(45) 


1  —  1/2 

where  a=  ( £?  +M)  E  .  The  quantity  under  the  radical  determines  the  nature 
of  loading.  There  are  three  possibilities: 

1  2  -  -  4*yC 

(1)  (f  C  +  M)  E  >  4<ryC,  i.e.,E> - 

<|C+M> 
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In  this  case,  the  ensuing  motion  is  decaying  without  oscillations. 


(2)  0  <  E  <  — g  ;  the  motion  is  decaying  with  oscillations. 

(3)  E  =  0.  This  case  was  shown  previously:  it  is  decaying  without 
oscillation. 

The  initial  conditions  for  the  solution  of  eo.uation  (45)  are: 

/  =  0  at  t  =  0 

(d¥)=io  +  ^1/2 

t=0 

where  uq  is  given  by  equation  (23). 

Thus  we  have  the  following  final  expressions: 

(1)  E> 

<ic*> 

4"  =5c°sh|t +Y  sinh  |t]e  *  J- 


and  5  ={|f -lSL  +  fC(E)1/2Y}-8inh  ft 

Hjf-  f  |  +  lC(E)l/2|}o08h  |  t 

-  C(E)l/2  (Ijf+l)  e2‘  (48) 


(2) 


0<  E < 


4FC 
,  2 
(|C+M> 


4=^ 


1 
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where 

*  =  </a2+4ayC  a=(|c  +  M)E1/2 

y=C(M®-ay)  Y  =  — ^-2- 

a 


*=  5yC  *  =  u0  +  icE1/2 

and  the  Integration  constants  A  and  6  are  given  by 


A  =  Jr  •  r~r 

P  COSO 

(49) 

and 

.  .  0  Y 

tano  =  ~y 

(50) 

C.  UNLOADING  WAVE  SOLUTION 

The  unloading  wave  is  fully  elastic,  unless  the  conditions  are  such  that 
secondary  yielding  in  the  opposite  direction  occurs.  The  assumption  is  made 
that,  immediately  after  the  sudden  unloading,  the  material  recovers  its  volume 

and  thus  Jjo  =  £  =  i.  The  procedure  of  formulation  and  solution  iu  cne  same  as 

for  the  loading.  Now,  however,  E  =  1  and  £  =  1  -  , 

To 

The  solution  is  given  by  the  following  equations: 

If  =  e  2  *  [j  cosh  1 1  +  Y  sinh  1 1]  -  %  (51) 

where,  for  Y,  £  =  -  1  -uQ. 
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f 


“ 3  frf-?  Y+lY]slnhfF 

+{l Y-  f  f + 

n.  i 

where  ip  =  a/oi  +4dy 
a  -  4  +M 
P  =  1 -<!  +  »! 
y  =  M  -  &! 


(52) 


(53) 


D.  CALCULATING  PROCEDURE 

The  velocity  equation  is  used  in  both  the  loading  and  unloading  cases,  to 
find  the  decay  of  the  velocity  and,  finally,  to  find  the  time  that  it  becomes  zero. 
This  time  is  then  substituted  in  the  strain  equation  to  find  the  resulting  strain. 


-r 

1 

I 
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SECTION  VI  -  EXAMPLE 


Consider  a  thick-walled  cylinder  with  ay  =  260, 000  psi  and  the  following 
dimensions  ( the  cylinder  is  filled  with  water): 


A  triangular  wave  of  the  following  form  is  transmitted  along  the  cylinder. 


This  wave  reflects  upon  arrival  on  a  closed  end.  Find  the  resulting  strains 

r  =  1  inch  and  r  =  3  Inches, 
o  o 

The  pressure  of  the  reflected  wave  front  is  approximately  twice  the  pressure 
of  the  incident  wave  front.  This  is  considered  the  highest  p--u»ible  internal  pres¬ 
sure;  thus,  the  calculations  will  be  performed  for  the  reflected  region  only. 

To  illustrate  the  wave  decay,  two  cases  are  analyzed;  Case  1,  where 
rQ=l  inch,  and  Case  2,  whore  rQ  =  3  inches.  Hence  P=-2(215,000)  =-430,000 
psi  at  rQ=  1  inch.  Using  equation  (23),  we  have 


US 


where 


u,  =  -  =  -  2864.6  in. /sec 

po  '■'o 

p0  =  0.29  lb/in.3 
CQ  =  200,000  in.  /sec 
g  =  32.2  ft/sec2 


From  equation  (34),  the  velocity  at  rQ  =  3  inches  (chosen  arbitrarily)  is 
U3  =  U1  /-JT  =-  2,864.6  yy=  -  1664.0  in./sec 


Thus  the  pressure  at  rQ  =  3  inches  is  given  by 


v  7Fx  Pi=-  248. 


280  pBi 


A.  CASE  1 


1. 


Loading  Wave 

From  equation  (28)  the  density  ratio  at  t  =  0  is  given  by 


O  C 

-  _S  o 

"  p  c0-u 


2.00x  10S 
2.00X106  -  2864.6 


=  .98568 


The  maximum  cr^  at  the  start  of  plastic  deformation  is  given  by 
°l=^  =  -^=-  357*l43PSi 


,uhere  ov  ~  -  250,  ono  psi 
M  =  0.7 

Thus  <  =  1  -  33?J  143  (1-2x0. 32)  =  1-  . 0119 x. 82  =  .99042 
30x10° 


at  ay  =  -  357  * 1-13  psi 

As  soon  as  the  material  begins  stretching  out,  the  loading  becomes  plastic. 
Considering  1  =  0  fron  equation  (38),  we  have 


45a 


w 

where  5^ 

Also,  Uj 

Substituting  in 


./dy-C  ■  71(10833  x .  98568  =  7.0082106  =  .09061 


250, 000 
30  x  106 


2864.6 
2  x  105 


=-.01432 


equation  (39),  we  have 


Since 


£  =  cosh  .0906  it  -  sin  it  .99061?  -  1 

u  =  ,  we  have 

u  =  .09061  (sinh  .  09061 1  -  .  15800  cosh  . 09061 t) 


First,  the  time  tQ,  where  0  =  0,  is  found  by  plotting  u  versus  t . 

(See  Figure  2).  The  £  versus  t  plot  is  shown  In  Figure  3:  it  1b  shown  that  as 
u  approaches  zero,  £  approaches  a  constant  value  of  0.0125.  The  radial  strain 

is  given  by  <  1  =  £  ~  -  1:  The  versus  t  plot  is  shown  in  Figure  4:  ^  ap¬ 
proaches  the  value  of  .  0264. 


2.  Unloading  Wave 

The  duration  of  the  triangular  wave  is  20  jjsec:  its  equivalent 
rectangular  profile  will  have  a  duration  of  10  jisec.  Hence, 

r  Co‘  2.0xl05xl0x  10'8  „  „ 

t  _  __  =  T  2.0 

It  is  shown  in  Figure  2,  however,  that  the  loading  process  terminals  at  t  =  1.75, 
Thus,  when  the  unloading  wave  occurs,  the  material  is  in  equilibrium. 

For  the  unloading  wave,  we  have  the  following  constants  (equation  (53))-. 
E  =  1 

C  =  1:  The  material  returns  to  its  initial  volume. 
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FIGURE  4 


f 


Hence  a  =  j  +  fi  =  0. 5  +  0. 303  =  .  803 

where  fi=  0.303.  N.B.  The  value  of  cr^  after  unloading  is  zero 
/}=  l~n-  0.697 
y=  0.303 
uQ=.  01432 

i  =•  -  1.0000  -  uo  =  -  1.0000  -  .01432  =  -  1.01432 

r 2 

I/)  =  7-*-  +4/3  =  1.71697 

|  =  .8585 

Thus  ^  =  ~^=  .43472  and  |  =  0.4015 

Hence  £  +  ®  (*)  =  -  1.01432  +  0.4015  ( .43472) 

=  -  1.01432  +.17454  =  -  .83978 

2 

Also,  -|x  |  =  .8585  x  .43472  =  .37321 

4x"4  =  .43472  x  .5  =  .21736 

1  fS 

|  Y  =  -  .4015  x  .978195  =  -  .392745 
^  x  |  =  . 43472  x . 4015  =  . 17454 

£=-  .978195  =  -  .489097 

2  —3 - 

and  |  jj  -  1  =  -  1  +  *21736  =  -  .78264 
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Hence  u  =  .  27686  sinh  .  8585 1  -  .  79696  cosh  .  8585 t  +  .  78264  e  * 4015 1 

and  g  =  e",4015t  {.43472  cosh  .8585t  -  1.01432  sinh  . 8585 1}-. 43472 

In  the  plot  of  u  versus  t  ,  Figure  5,  the  time  when  the  velocity  becomes  zero 
is  shown  as  t  =  0.026  (t=  .26x10'®  sec).  Substituting  in  the  equation  of  g  we 
find  the  residual  hoop  strain,  £  -  .01713,  Equations  (51)  and  (52)  were  derived, 
however,  with^0  =  0  at  t  =  0  (with  reference  to  the  equilibrium  conditions  after 
the  loading). 

To  find  the  residua1,  hoop  strain  after  unloading,  we  have: 

^  ^  ~  ^"unloading  ~  ^loading 

=  .01713  -  .0125  =  .  00463  (compression) 

In  order  to  find  the  residual  hoop  stress  we  use  the  principle  of 
elastic  recovery.  Thus,  the  volumetric  strain  ey  is  given  by: 

ev-cl  +  ^unloading 

where  is  the  radial  linear  strain  change  upon  relief  and  (g  is  the  corresponding 

hoop  strain.  We  know,  however,  the  radial  unloading  stress:  i.e.,  Aa1  =  430,000 

psi.  Hence,  «.  =  -30> 3  =  .0143. 

1  30x10° 

The  value  of  ey  is  also  known  from  the  loading  conditions: 

-ev=  ¥-=~  f  +1  =  +  -°1813 

o  Ko 

Thus,  the  hoop  strain  for  zero  stress  (permanent  set)  is  given  by 
Q  =  ey  -  ex  =  .01813  -  ,0143  =  .C036. 

The  corresponding  permanent  set  in  the  radial  direction  is 
tQ  = .  026  -  .  0143  =  .012 .  In  the  calculations  for  the  unloading  wave,  it  was 
assumed  that  upon  relief  the  recovery  was  instantaneous:  thus  —  =  1  at  t=0. 


51 
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FIGURE  5 


Then,  for  t  >  0,  -p—  decreases,  clianging  born  hoop  and  radtal  strains  according  to 
o 

to  the  following  equations: 

S  =  l  -  jr  and  €x  =  C  ^  -  1. 

Hence,  the  material  will  be  under  the  following  strain  conditions  : 


Without  pursuing  this  case  any  further,  it  is  easily  concluded  that  secondary 
yielding  will  take  place.  This  is  undesirable  for  most  designs. 


B.  CASE  2 


Ifla&aaL&ayfi. 

The  velocity  at  r0  =  3  inches  was  given  as: 

u,  =  -  1654.0  and  u„  =  Q  =  -  -  *00827 

“3  3  o  2x10^  *— 

^  u 

C  =  “7 —  =  1  -  =  .99173,  and  a  =  250*000  =  36?)  l43  j 

Co  0  ^(elastic)  0,7 


Hence,  w=  00833  x  .99173  -  ^.0082611  =  .0908906 
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.  \=_  .00827  =  -  .0909885 

”  u  '  .0908906 

Hence  /  =  cosh  . 0908906 t  -  .0909885  sinh  . 0908906 t  -  1 
and  u  =  .0908906  {sinh  .  0908906 1  -  .0909885  cosh  .0908906  } 

The  duration  of  the  wave,  on  the  other  hand,  is  given  by 

r  _  10  X  10~6  x  2  x  135  _  2 
3  3 

where 

t  =  10  x  10"6 

C  =  2  x  10® 

o 

r  =3  inches 
o 

The  u  versus  t  and  £  versus  t  plots  are  shown  in  Figures  6  and  7, 
respectively.  Another  case  of  ay  =  150, 000  psi  is  also  shown  for  comparison. 

With  the  value  of  t  at  rQ  =  3  inches  being  smaller  than  at  rQ  =  1  inch, 
it  is  shown  that  the  unloading  wave  occurs  before  the  material  achieves  equilibrium. 

2.  Unloading  Wave 

Following  the  same  procedure  as  for  rQ  =  3  inches,  with 
u0  =  0.00827  -  .00275*  =  .00552 

£  =  e’,4015t  {.43472  cosh  ,8585t  -  .96794  sinh  .85851  )  -  .43472 

-  -  401 S  t 

anH  u  =  -  .  27793  sinh  .  8585 1  +  .  78816  cosh  .  8585 1  -  .  78264  c* 

The  plot  of  Q  versus  t  is  shown  in  Figure  8.  Comparing  this  un¬ 
loading  with  that  of  rQ  =  1  inch,  we  see  that,  at  rQ  =  3  inches,  equilibrium  is 
achieved  in  about  half  the  time. 


*  velocity  from  the  loading  process 
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The  corresponding  value  of  £  at  t  =  .01  is  £  =  .00727.  The  £  value 
at  t  =  0.667  for  loading  is  £  =  .00365.  Hence,  the  resulting  strain  is  (.00733 
-  .00365)  or  £  =  .00368  (compression). 


To  find  the  resulting  stresses,  the  permanent  set  in  both  the  hoop 
and  radial  directions  must  be  found.  Employing  the  concept  of  elastic  recovery 
again,  we  have 

ev=  *•*«! 


( 


1 


248.280 

30x10® 


. 005276 


Furthermore,  -  e  =  +  +  1  =  -  .99173  +  1  =.00827. 

vo  "o 

Thus,  £ Q  =  0:  This  result  indicates  that  there  is  no  permanent  set  in  the  hoop 

direction.  In  the  radial  direction,  =  1  +  /  =  1+  .00366  =  1.00365  from 

ro 

r  P 

loading.  Hence  ~  =  j  Q0^6R  =.99636.  Also,  ~  =  .99173. 


=  .99636  x  .99173  -  1  =  .98812  -  1  =  .01188  from  loading. 


During  the  unloading,  there  is  the  elastic  recovery  where  <1=  .00827: 

then,  there  is  the  further  change  of  e,  due  to  the  change. 

1  ro 

r 

Hence  c,  =•—  -  1  +  .00827 
l  r 

—  =1-^  =  1  -  .00727  =  .99373 
r 

=  1.0063 
r 

and  ex  =  .  0063  +  .  00834  =  .  01464 
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Thus,  change  in  f  i  =  .01464  -  .0119  =  ,00274  (tension).  The  value 
of  Cj  when  the  radial  stress  is  zero  is  .01203  -  .00827  =  .00376, 

Residual  radial  stress  =  .00376  +  .0026  =  .00636  or 
.00636  x  3  x  106  =  190.800  psi  (tension). 

g 

The  residual  hoop  stress  =  .00368  x  3  x  10  =  110,400  psi 
(compression). 


f  .00368 

1  1  j 

r° 

‘  i  i  i 

-  “4  / 

+.005 

/ 

100,000 

It  .  . 

tio,  OOO  PSI 

- 

200,000 

or  191,000  +  110,000  =  301,000  psi.  Since  cr^  was  assumed  at  250,000  psi, 
the  material  will  go  through  secondary  yielding. 
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SECTION  VII  -  DISCUSSION 


The  theory  of  unsteady  waves  was  used  to  develop  working  equations  for 
computing  the  radial  and  hoop  strains  in  hollow  cylinders  under  transient  internal 
loads.  The  wa  *  propagates  with  a  steep  elastic  front  which  is  followed  by  an 
attenuated  plastic  wave.  The  changes  of  density  and  velocity  immediately  behind 
the  front  are  instantaneous:  furthermore,  the  plastic  wave  is  not  attenuated  at 
t=0.  From  this  pattern  it  is  easily  seen  why  the  stress  field  of  the  material  at 
t=0  is  three-dimensional  and  the  shear  stress  considerably  higher  than  for 
slowly  applied  loads.  Under  hydrostatic  pressure,  at  t=0,  the  hoop  stress  is 
compressive  (for  a  compressive  wave).  As  soon  as  the  deceleration  process 
starts,  however  (t  >0 ),  the  material  is  under  tension  in  the  hoop  direction. 

The  change  from  compressive  to  tensile  hoop  stress  does  not,  however,  take 
place  instantaneously.  Thus,  the  Btress  behind  the  wave  remains  elastically 
higher  than  when  the  load  is  applied  slowly  until  the  transition  from  compressive 
to  tensile  stress  is  made.  This  may  be  the  reason  why  higher  yield  stresses 
are  manifested  during  impulsive  loads. 

Another  important  aspect  of  the  analysis  is  the  duration  of  the  wave  as  it 
propagates  through  the  wall  of  the  cylinder.  The  dimensionless  parameter 
tc 

t  = -  is  constant  for  all  r  .  Thus,  as  r  increases,  the  time  t  to  achieve 

r  oo 

o 

steady  state  also  iccrc s.y-e.  The  results  are  shown  in  the  example  where  the 
cuoco  ior  rQ  -  1  and  rQ  -  3  are  analyzed.  At  rQ  =  1,  the  material  wi'l  ex¬ 
perience  plastic  flow  first  in  the  outward  direction  and  then  in  the  inward  direc¬ 
tion  (secondary  yielding):  at  rQ  =  3,  the  same  events  take  place  but  to  a  lesser 
degree.  If,  however,  the  cTy  had  the  value  of  290, 000  psi  the  material  could 
withstand  the  pressure  at  r0  =  3  without  secondary  yielding.  Thus,  alter  the 
first  pulse  the  material  will  withstand  the  loads  for  all  subsequent  pulses  without 
further  yielding. 
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In  his  tests  with  thick-walled  cylinders,  Rinehart^  has  shown  that  the 
mode  of  failure  near  the  bore  was  in  shear  while  farther  out  from  the  bore  the 
material  fractured  in  tensile  hoop  stress.  This  can  be  attributed  mainly  to 
the  high  hydrostatic  pressure  at  t  =  0:  as  the  wave  amplitude  decreases,  the 
hydrostatic  pressure  decreases,  allowing  the  material  to  be  predominantly  in 
hoop  tension. 
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CONCLUSIONS 


1 .  Extrusion  of  an  appropriate  billet  from  an  ultra  high  hydrostatic 
fluid  pressure  medium  (preferably  inhibited  water)  is  theoretically 
possible  from  the  effects  of  a  superposition  of  a  detonation  shock  type 
pressure  pulse  created  from  the  rapid  discharge  of  capacitor  stored 
energy  across  a  suitable  gap. 

2.  As  the  relationship  uo"  =■  ^ In derived  from  the  variables 

Po  t 

associated  with  electrohydrnulic  extrusion  bears  a  striking  similarity 

to  the  well  known  empirical  expression  for  extrusion,  P  =  o  In  4s1 

A 

where  P  =  pressure  to  extrude,  a  =  resistance  to  deformation  and 
An 

In  "a —  *  log  of  the  extrusion  ratio,  a  logical  analytic  approach  to 
the  proposed  extrusion  experiments  is  available. 

3.  Extrusion  for  one  pulse  is  dependent  upon  the  shock  pressure 
pulse  amplitude.  The  final  extrusion  length  will  be  dependent  upon 
the  rate  of  pulses  generated  from  a  series  of  capacitor  discharge 
events, 

4.  A  pressure  vessel  can  be  designed  to  contain  a  200,  000  psi 
hydrostatic  fluid  pressure  with  a  superposed  dynamic  detonation 
pressure  pulse  loading  of  425,  000  psi  after  reflection. 

5.  At  hydrostati''  pressures  of  200,  000  psi,  the  particle  velocity  in 
the  walls  of  the  container  as  a  consequence  of  the  capacitor  discharge 
event  attenuates  rasidly  from  a  critical  value  to  a  value,  which  con- 
siacieci  with  the  static  stress  pattern  is  within  the  elastic  prope-*y 
limits  of  accepted  container  design  materials. 

6.  From  ihe  pressurized  fluid  boundary  outward  into  the  container 
walls  for  a  distance  of  2  inches,  the  behavior  of  the  metal  is  critical 
because  of  high  particle  velocity  and  cyclic  loading  conditions.  For 
this  reason,  the  pressure  vessel  design  must  include  a  replaceable 
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liner  of  ultra  high  strength  material.  This  liner  must  be  put  into  a 
condition  of  positive  axial  and  radial  restraint. 
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FUTURE  WORK 


A  purchase  order  based  on  analysis  of  vendor  proposals  will 
be  issued  by  Republic  Aviation  Corporation  for  the  procurement  of 
the  pressure  vessel  and  pumping  system  best  suited  to  the  objectives 
of  the  program. 

The  derived  equations  involving  the  radial  propagation  of 
stress  waves  in  thick-walled  cylinders  under  transient  detonation 
internal  loads  along  with  conventional  static  stress  criteria,  and 
extrusion  experimental  requirements  will  be  applied  to  the  electro- 
hydraulic  pressure  vessel  design  problem  to  formulate  a  final  optimum 
working  design.  From  these  engineering  considerations,  assembly 
drawings  will  be  produced  by  the  selected  vendor. 

Upon  approval  by  Republic  Aviation,  the  assembly  drawings 
will  be  worked  into  detailed  design  drawings  suitable  for  manufacture 
of  the  equipment.  Fabrication  of  the  pumping  equipment  will  commence. 

An  investigation  into  presently  attainable  charging  rates  and 
methods  of  rapid  cyclical  release  of  high  voltage  capacitor  discharge 
energy  will  be  made  to  obtain  a  practical,  economical  approach  to 
the  equipment  best  suited  to  the  requirements  of  clectrobydraulic 
extrusion. 
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